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Exchange striction induced giant ferroelectric polarization in copper based 

multiferroic material Q!-Cu2V207. 
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We report Q-CU 2 V 2 O 7 to be an improper multiferroic with the simultaneous development of 
electric polarization and magnetization below Tc = 35 K. The observed spontaneous polarization 
of magnitude 0.55 is highest among the copper based improper multiferroic materials. Our 

study demonstrates sizable amount of magneto-electric coupling below Tc even with a low magnetic 
field. The theoretical calculations based on density functional theory (DFT) indicate magnetism 
in q:-Cu 2 V 207 is a consequence of ferro-orbital ordering driven by polar lattice distortion due to 
the unique pyramidal (CuOs) environment of Cu. The spin orbit coupling (SOC) further stabilize 
orbital ordering and is crucial for magnetism. The calculations indicate that the origin of the giant 
ferroelectric polarization is primarily due to the symmetric exchange-striction mechanism and is 
corroborated by temperature dependent X-ray studies. 

PACS numbers: 75.85.-|-t, 71.20.-b 


Recently multiferroic materials with mutually coupled 
ferroelectric (FE) and magnetic orders have attracted 
considerable interest for their versatile technological 
as well as fundamental importance. A strong 

magneto-electric (ME) coupling is expected in improper 
magnetic mutiferroics where ferroelectricity is induced 
by a specific magnetic order. In the last one decade, 
several magnetic multiferroics have been discovered 
[113 where FE polarization is either associated with 
spiral magnetic structure induced by spin-orbit coupling 
(SOC) [111, ll2| or by symmetric exchange striction (SES) 
mechanism in case of collinear magnets. 0, Hi Due to 
the secondary nature of the electric order, the value 
of the FE polarization in such magnetic multiferroics 
is much smaller (generally 0.01 ^C.cm“^) compared 
to the ‘proper’ FE. 01 A recent breakthrough in this 
direction is the discovery of giant ferroelectricity 
0.3 ^C.cm“^) and large ME coupling in mixed valent 
manganate CaMn70i2 below about 90 K mediated 
by both Dzyaloshinski-Moriya (DM ) in teraction as well 
as exchange striction mechanism. [l5| In this respect 
cuprates may be an attractive option as the orbital 
degrees of freedom and strong Coulomb correlations 
present in cuprates may not only lead to lattice dis¬ 
tortion and magnetism but also possibly induce a 
coupling between them which are essential ingredients 
for multiferroicity. 


In view of the above, we investigated the Cu-based 
oxide CU2V2O7 in its orthorhombic a phase. CU2V2O7 
crystallizes in at least three different polymorphs, namely 
a, j3 and 7-phases where only the a phase is non- 
centrosymmetric [IM3 and is important in the present 
context. It consists of magnetic Cu^+ (3d®, S = ^) and 
nonmagnetic (3d ®, S = 0) metal ions making it a sys¬ 
tem having both partially filled and empty d shells sim¬ 
ilar to BiFeOa, BiMnOg, Pb(Fe2/3Wi/3)03 etc. [HllOl. 


All Cu2+ ions are equivalent with fivefold coordination 
to oxygen atoms forming a distorted [CuOs] polyhedron. 
Each Cu-polyhedron is linked with another two via edge 
sharing and they together form two sets of mutually per¬ 
pendicular zig zag chains (see Fig. 1). These chains are 
separated by Y 20 ^~ anionic group resulting from the two 
corner sharing VO4 tetrahedra. [ibl l2Ij 

The magnetic behavior of a-Cu2V2 07 have been in¬ 
vestigated earlier on polycrystalline samples. It 

was reported that a-Cu2V207 is an antiferromagnet with 
weak ferromagnetism at low temperature. A magnetic or¬ 
der was seen below 35 K (Tc) accompanied by a change 
in slope in dielectric response near Tc- In this letter, a 
combined theoretical and experimental work establishes 
that the compound is an improper multiferroic with giant 
P and ME effect where the origin of giant FE polarization 
is primarily due to symmetric exchange-striction mech¬ 
anism, and the magnetism is stabilized by ferro-orbital 
ordering. 

The experimental studies including magnetic, dielec¬ 
tric and electric polarization measurements were per¬ 
formed on a well characterized sintered polycrystalline 
sample, [20 which have been described in detail in Sup¬ 
plemental Material (SM). All the electronic structure 
calculations presented in this paper are performed us¬ 
ing DFT within local density approximation (EDA) and 
projector augmented wave (PAW) method as encoded in 
Vienna ab-initio simulation package (VASP) [25I - E0 (see 
SM). 

Fig. 2 (a) describes magnetization (M) vs. tempera¬ 
ture (T) data in zero-field-cooled-heating (ZFCH), field¬ 
cooling (FC) and field-cooled-heating (FCH) protocols 
under magnetic field H = 100 Oe. M(T) shows a sharp 
rise at Tc = 35 K indicating the transition to a magnet¬ 
ically ordered state. The thermal hysteresis between FC 
and FCH around Tq indicates the first order nature of 
this transition. The inverse molar susceptibility (x“^(T), 
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FIG. 1. (a) Cu atoms in the conventional unit cell form a 
pair of mutually perpendicular zig-zag chains. Various spin 
exchange interactions in a-Cu 2 V 207 are also shown, (b) Edge 
sharing Cu-polyhedron forming two mutually perpendicular 
chains connected by two corner sharing VO 4 terahedra. (c) 
Change in the nearest neighbor Cu-0(4)-Cu and Cu-0(3)-Cu 
pathways as well as the next nearest neighbor exchange path 
upon relaxation, (d) Three dimensional electron density plots 
showing the ferro-orbital order within LSDA+U-l-SOC. 


where x = M/H) obeys Curie-Weiss law above 80 K 
(see inset of fig.2(a)) and we get Curie-Weiss temperature 
{6c) to be « — 78 K which indicates the predominant an¬ 
tiferromagnetic (AFM) correlations in the system. The 
effective moment of Cu^’*' is « 1.92 and it is slightly 
higher than the spin-only value (= 1.73 ^b)- At low T, 
M almost saturates which is not a likely behavior of a 
pure AFM ordering. Possibly, the magnetic state below 
Tc is canted AFM type. The isothermal M vs. 77 at 5 
K for 77 = ± 9 kOe is shown in the main panel of Fig. 2 
(b). The curve shows clear hysteresis which reaffirms the 
presence of ferromagnetic (FM) component. The coerciv- 
ity of the loop is found to be about 2 kOe. The full loops 
both at 5 K and at 150 K (well above Tc) are shown in 
the inset. The M — H curve at 5 K, however, does not 
show full saturation even at 50 kOe, and it once again 
indicates canted spin structure. We can ht the high held 
data (between 77 = 30-50 kOe) with an empirical rela¬ 
tion M{H) = XafmH -b Ms, where XafmH is the linear 
term due to AFM component and Ms is the saturation 
magnetization due to the FM part. We Hnd Ms to be 
0.08 //s/f.u., which is quite small compared to the full 
saturation moment of two Cu^"*" (^ 2 hb) ions indicat¬ 
ing the presence of weak ferromagnetism. Interestingly, 
the 5 K isotherm is not found to be quite smooth, and 
it contains signature of sharp jump whenever the held 
changes its sign. This may indicate the presence of uni¬ 
axial anisotropy in the system. 

Fig. 3 (a) shows the T-variation of the real part of the 
complex dielectric permittivity (e') measured at different 
frequencies (/). e' is almost constant and independent of 



FIG. 2. (a) The T dependence of ZFGH, FC and FGH mag¬ 
netization data of Q-GU 2 V 2 O 7 . The inset shows 1/x versus 
T along with the fitting of Curie-Weiss law. (b) indicates the 
isothermal M vs H data at 5 K and 150 K. The main panel of 
(c) shows the T variation of orthorhombic lattice parameter b 
along with inset depicting the T dependence of lattice volume 
V obtained from powder XRD data, (d) shows the change in 
bond lengths di and 62 corresponding to the interactions Ji 
and J 2 (see fig. 1 (a)) with T. The inset shows change in the 
GuOs polyhedron after relaxation as obtained from theory. 


/ in the low-T regime (below ~ 130 K) which signifies the 
static dielectric constant originating from the intrinsic 
contribution [^. A closer look at low-T part shows the 
existence of a small but clear hump-like anomaly around 
35 K coinciding with Tjv. It is free from any frequency 
dispersion, suggesting that this feature may be related 
to some long range electric order. The imaginary part 
of permittivity is quite small (particularly below about 
130 K) indicating that the sample is highly resistive and 
the estimated resistivity at 50 K is found to be ^ 100 
Mfl-cm. 

Considering the electric anomaly near the magnetic 
transition, it is tempting to measure the magneto¬ 
dielectric properties of the sample [s^. Fig. 3 (b) shows 
the T variation of e' measured at 77 = 0 and 9 kOe. 
Clearly, e' shows significant effect of magnetic field below 
about 80 K. In the inset of fig. 3 (b), we have plotted 
the change in e' as a function of T due to the applica¬ 
tion 9 kOe of field and we observe a significant value of 
magneto-dielectric effect (as high as 3.5%) at around 30 
K. This is remarkably large considering the small value 
of the applied field . 

In order to shed more light on the nature of the hump¬ 
like feature observed in the e'{T) coinciding with the 
magnetic anomaly at Tc, we measured pyroelectric cur¬ 
rent {Ip) (see SM) after cooling the sample from room 
temperature with different electric field {Ecooi). From 
the T variation of Ip, spontaneous polarization has been 
calculated (see fig 3 (c)), which we denote by P/. Clearly 
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Pj shows a sharp increment below about 35-40 K, even¬ 
tually saturating at a lower T. The magnitude of Pi 
clearly increases with the cooling field. We also mea¬ 
sured Pi with different polarity of Ecooi and Pi changes 
sign depending on the chosen sign of Ecooi (see inset 
of fig. 3 (c)). Such behaviors of Pi confirm that the 
sample undergoes long range FE order below 35 K with 
the development of spontaneous polarization. Since, the 
electric order is concomitant with the magnetic order, 
the sample can be assigned as an improper multiferroic 
material. It is to be noted that even at room tempera¬ 
ture a-Cu 2 V 207 possesses a non-centrosymmetric crys¬ 
tal structure with polar point group (mm 2 ), which in 
general belongs to the pyroelectric class of crystals. [31 1 
However, a switchable spontaneous Pi is only achieved 
below Tc, possibly arising from the favorable lattice dis¬ 
tortion associated with the magnetic order. Remarkably, 
from the pyroelectric measurement the saturation value 
of the spontaneous Pi is found to be as large as 0.55 
/rC.cm“^, which is substantially high compared to the 
other copper based magnetically driven ferroelectrics to 
date, [ii,!!! 

The ferroelectricity is further confirmed by the mea¬ 
surement of electric polarization versus electric field 
{P — E) loop using a FE loop tracer as shown in Fig. 
3 (d). The data recorded above Tc (50 K and 80 K) 
do not show any loop, whilst clear loops with tendency 
for saturation are observed at 10 K and 30 K (which are 
below Tc). Observation of such prototypical hysteresis 
loop in P is an essential proof for the development of FE 
state below Tq. At the value of E = 1 kV.cm“^ hystere¬ 
sis almost closes and we get a value of polarization close 
to 1 ^C.cm“^. 

The experimental results discussed above lend support 
to the fact that ferroelectricity is induced by magnetism 
in a-Cu 2 V 207 resulting in substantially large magneto¬ 
dielectric response. It is also interesting to note that 
the magnetic transition at Tq is first order in nature 
which indicates possible structural transition associated 
with the magnetic as well as electric orderings. In view 
of the above, both exchange striction as well as inverse 
DM effect may be the likely mechanism for the giant FE 
polarization in a-Cu 2 V 207 . The DM interaction is not 
unlikely here as the obtained effective moment per Cu 
site (1.92 fis) is bit higher than the spin only moment 
(1.73 ^b) presumably due to the orbital contribution of 
the magnetic moment as a consequence of finite SOC. It 
is to be noted that the spontaneous polarization in a- 
CU 2 V 2 O 7 is about one order of magnitude higher than 
the other spiral DM type multiferroics (such as TbMnOs, 
where P ^ 0.05/rCcm“^). Therefore, the exact origin of 
FE state in this compound may involve more complex 
mechanism. 

To understand the electronic structure, the exchange 
mechanism and the origin of ferroelctric polarization, 
we have performed first principles DFT calculations, 


using VASP as described in SM. The non-spin polarized 
electronic structure calculations reveal that the oxygen 
p-states are completely occupied while the vanadium-d 
states are empty and the Fermi level is hosted by half- 
filled predominantly Cu-dx 2 -y 2 states consistent with 
the Cu 2 ''’V 2 ~'' 07 “ nominal ionic formula for the system. 
The distortion of the CuOs polyhedra is triggered by the 
orbitally active Cu^“'" ion in such a way that Cu-dx 2 -y 2 
states are well separated from the rest of the Cu-d 
states thereby promoting ferro-orbital ordering. Such 
a ferro-orbital ordering will favor anti-ferromagnetic 
coupling between the Cu ions in the chain. 3^ It is 
also interesting to note that the distortion is such that 
the nearest neighbor (NN) Cu ions in each chain are 
now coupled by two asymmetric bonds Cu-0(3)-Cu 
and Cu-0(4)-Cu forming a non-centrosymmetric Cu 02 
plaquette (see Fig. 1(c)) which will not only favor DM 
interaction but also stabilize ferroelectric polarization. 


To account for the observed magnetism in this system, 
we have evaluated various symmetric spin exchange 
interactions between the Cu-ions by performing to¬ 
tal energy calculations in the framework of LDA -|- U 
method, (see SM) The value ofUe/f {U — J) was taken to 
be 6.5 eV following the usual choice for the cuprates, (ssj 
Constraining the range of interaction to 5.42 A, we 
calculated five dominant exchange interactions (the 
various spin-exchange paths are shown in Fig. 1(a)). 
The dominant inter-chain exchange interaction J 3 (- 
13.61 meV) is antiferromagnetic, followed by intra-chain 
Ji(-4.67 meV) and interchain J 2 (4.07 meV) which are 
AFM and FM respectively. Other exchange interactions 
J4(0.26 meV) and J 5 (2.37 meV) are small and FM. The 
AFM exchange interactions Ji and J 3 are mediated via 
Cu-O-Cu and Cu-O-V-O-Cu paths respectively [See Fig. 
1(b)]. The larger bond angles in the J 3 exchange path 
make this interaction stronger compared to Ji. The 
NN AFM Ji is consistent with the ferro-orbital order 
and the strong inter-chain coupling J 3 identified by our 
first principles calculations is responsible for the long 
range magnetic order. The FM exchange interaction J 2 
is mediated by the exchange path Cu-O-O-Cu where 
the two O atoms that mediate the exchange interaction 
between Cu ions are at an angle (ZCu-0-0) of 70.07° 
thereby favoring FM interaction. Using the computed 
exchange parameters, the Curie-Wiess temperature 9c 
for a-Cu 2 V 207 is calculated to be -77.4 K which is 
remarkably close to the experimental value (-78 K). 

In the absence of any spin frustration, the canted 
AFM suggested by the M — H curve may be attributed 
to the DM-type interaction due to SOC. We have 
therefore considered the antisymmetric part of the spin 
Hamiltonian H = ^ ^j) calculated the 

three components D^, of the DM parameter up 

to 3rd nearest neighbor for a-Cu 2 V 207 by performing 
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FIG. 3. (a) shows the T variation of the real part of dielec¬ 
tric permittivity, (b) shows the e data measured at = 0 
and 9 kOe along with the inset showing the change in e due 
to H. (c) represents the temperature dependence of sponta¬ 
neous electric polarization, Pi calculated from the pyroelec¬ 
tric current measurements. Insets of the (c) represents P 
with positive and negative cooling electric fields, (d) shows 
the polarization hysteresis (P — E) loop measured at different 
temperatures. 


LDA -1-17-1- SOC calculations. The ratio |^| ~0.5 
suggests unusually strong DM interaction in a-Cu 2 V 207 
very similar to that reported for CaMn 70 i 2 system 
[ISj. The calculation of magnetocrystalline anisotropy 
energies reveal 6-axis of the conventional unit cell to be 
the easy axis. 


Guided by the exchange interactions and in view 
of the importance of the SOC, we have considered 
three magnetic configurations namely FM, AFM and 
non-collinear (NONC) as shown in Fig. 6 of SM. In 
our model NONC magnetic configuration, the nearest 
neighbor spins in each zig-zag chain are antiparallel in 
the b-direction with a very small component along a and 
c direction as a result of canting. Our calculations (see 
Table-II of SM) reveal SOC stablizes magnetism and the 
non-collinear magnetic structure is energetically degen¬ 
erate with a very small magnetic moment(0.17//B/f.u.). 
A three dimensional plot of the spin density for the 
NONC structure with LDA-I-U-I-SOC conhrm ferro 
orbital ordering (see Fig. I (d)) where the same orbital 
is occupied at each site but rotated with respect to 
each other due to the distorted CuOs polyhedron. 37 


Finally the FE polarization is calculated for the above 
mentioned three magnetic configurations using the Berry 
phase method [s^ as implemented in VASP. 271 For the 
purpose of understanding the contribntion of SOC and 
exchange-striction on ferroelectric polarization, calcula- 


TABLE I. Calculated Polarization in various magnetic struc¬ 
tures 


Config. 

AE 

(meV) 

Polarization due to 

SOC 

APsoc 

(pC.cm“^) 

Polarization due to 

Exchange striction 
APea: 

(pC.cm“^) 

FM 

29 

0.02 

3.72 

AFM 

0 

0.05 

4.08 

NONC 

0 

0.01 

4.03 


tions are carried out including SOC both for the experi¬ 
mental structure and the relaxed structure. The results 
of our calculation (Table I) although suggest the impor¬ 
tance of SOC but conclusively establish that exchange- 
striction is the primary mechanism for the giant ferro¬ 
electric polarization for this system. As shown in Fig. 
2 (d) (inset), as a result of exchange striction there is a 
compression of the CuOs polyhedron, resulting in shorter 
bond length in the exchange path Ji and consequently 
increase in the bond length in the exchange path J 2 (See 
Table II of SM). Shorter bond length di in the AFM ex¬ 
change path Ji will enhance the hopping contribution C 
and add to the gain in energy by the superexchange (oc 

tyu). 

In order to corroborate our theoretical results on 
exchange-striction we performed T-dependent X-ray 
diffraction (XRD) measurement on the sample . Our 
analysis of the data indicates that the crystal symme¬ 
try (Fdd2) remains unchanged both below and above 
Tc- l4Q| However, there is clear change in the orthorhom¬ 
bic lattice parameters (a, 6 and c) at Tc, where the 
change is prominent for lattice parameter b (fig. 2 (c)). 
Remarkably on cooling across T/v there is a sharp first 
order like change in di and c ?2 (see hg. 2 (d)) with Ac?i 
= -0.085(2) Aand Ad 2 = 0.159(2) A, while d^ almost re¬ 
mains unchanged (Ada =-0.01 A) in excellent agreement 
with our theoretical prediction. (See Table II of SM) 

In conclusion, we have found that a-Cu 2 V 207 is a 
magnetic multiferroic material with the highest FE po¬ 
larization among the known Cu based multiferroic ox¬ 
ides with sizable amount of ME coupling. q;-Cu 2 V 207 
turns out to be an unique example of multiferroic ma¬ 
terial with single valent Cu^'*' ions where orbital degrees 
of freedom lead to polar distortion and ferro-orbital or¬ 
dering favoring AFM. The SOC further assists to stabi¬ 
lize orbital ordering and magnetism. Finally the AFM 
interaction promotes large exchange-striction and is the 
primary mechanism that gives rise to giant electric po¬ 
larization. 
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Supplementary Materials (SM) for 

Exchange striction induced giant ferroelectric 
polarization in copper based multiferroic 
material a-Cu2V207 


EXPERIMENTAL ASPECTS 


AC dielectic measurements 

The ac dielectric measurements were performed using 
an Agilent E4980A precision LCR meter in the tempera¬ 
ture range 5-300 K in a helium closed cycle refrigerator. 
An electromagnet with maximum field strength of 9 kOe 
(at 5 cm pole separation) was used to apply magnetic 
field during dielectric measurements. 


Sample Preparation and Characterization 


Pyroelectric Measurement 


Polycrystalline sample of CU 2 V 2 O 7 was prepared by 
conventional solid state reaction route in air. Highly 
pure CuO and V 2 O 5 were mixed thoroughly in a sto¬ 
ichiometric ratio and homogenized with ethanol in an 
agate mortar. The mixture was pressed into pellets and 
sintered at 600° C for 90 h with several intermediate 
grindings. We found that any high temperature sinter¬ 
ing results in the formation of the /3 phase of CU 2 V 2 O 7 . 
Powder X-ray diffraction (XRD) patterns at room tem¬ 
perature (300 K) as well as at 200 K, 100 K, 70 K, 40 
K, 25 K and 15 K were recorded using Cu Ka radiation 
in Rigaku-TTRAX-III diffractometer which ensure the 
single phase of the polycrystalline sample with no de¬ 
tectable secondary phase. Reitveld refinement analysis 
were performed with the help of MAUD software on the 
powder XRD patterns. The XRD patterns at all tem¬ 
peratures (both above and below Tq = 35 K) are well 
fitted in orthorhombic crystal structure with Fdd2 space 
group where reliability factors {a) are < 1.4. The lattice 
parameters at 300 K were found to be a = 20.692 k,h = 
8.413 Aand c = 6.451 A. The Cu-O-Cu and V-O-V angles 
are found to be 106° and 148° respectively. These angles 
as well as Cu-0 and V-0 bond lengths agree well with 
the previous work on the structure of Qf-Cu 2 V 207 . [l| 
Although no change in crystal symmetry was observed 
(within the resolution of our data), clear change in lat¬ 
tice parameters and bond lengths are present around the 
Curie point. The deviations in the lattice parameters and 
in the bond lengths with temperature are well outside the 
error limit of refinement which are respectively « 0.0006 

Aand Ri 0.002 A. 


Magnetic measurements 


The pyroelectric current of the sample was measured 
using a Keithley Electrometer (model 6517B) in the he¬ 
lium closed cycle refrigerator. From the pyroelectric cur¬ 
rent polarization have been calculated (see fig. 4). In 
order to record Ip, a capacitor type arrangement was 
used where a pair of electrodes were attached to two flat 
surfaces of the pelletized sample using silver epoxy. The 
sample was first cooled down to 10 K in presence of elec¬ 
tric field {Ecooi = 2, 3, 4 and 5 kV/cm). After reach¬ 
ing lowest T, Ecooi was set to zero and Ip was mea¬ 
sured duriM heating of the sample at a constant rate of 
4 K/min. The pyroelectric current density(Jp) was 
calculated by dividing /p by the area of the electrode(A). 
We can calculate P with the following relation: 


Pi 




IpdT 


, where ^ is the rate of change of temperature. Here we 
assume that Pj vanishes above 45 K. Clearly a peak is 
noticed around 35 K and the value of Jp increases with 
the increase of Ecooi- The peak position is almost cool¬ 
ing field independent which indicates that the intrinsic 
pyroelectric current is dominant. 


Polarization hysteresis loop measurement 

Polarization hysteresis loops (P — E) at different 
constant temperatures were measured by a Ferroelec¬ 
tric Loop Tracer from Radiant Technology (Precision 
Preinier-H) in a helium closed cycle refrigerator. Mea¬ 
surements were performed in presence of electric field as 
high as ± 1 kV/cm with a time period of 10 ms in stan¬ 
dard bipolar mode. 


The magnetic measurements were performed on a 
Quantum Design SQUID magnetometer (MPMS-4, Ev- 
ercool). Extra precautions were taken to judge the small 
thermal hysteresis observed in the isothermal magneti¬ 
zation measured at 5 K where data were recorded with 
stable mode option. 


THEORETICAL ASPECTS 
General techniques 

All the electronic structure calculations presented in 
this paper are performed using first principles density 




FIG. 4 . T dependent pyroelectric current density of alpha- 
CU2V2O7 measured after cooling the sample in different elec¬ 
tric fields. 



FIG. 5. Effective Cu-d^ 2 _j ,2 Wannier function plot showing 
the exchange path for the inter-chain exchange interaction J 3 . 

functional theory (DFT) within local density approxi¬ 
mation (LDA) including Hubbard U 14 using projector 
augmented-wave (PAW) method |g encoded in the 
Vienna ab initio simulation package (VASP) [ 3 , 0 |- The 
values of on-site Coulomb interaction (U) and Hunds 
rule coupling (Jp) parameters were taken as f/=7.5 
eV, J//=1.0 eV [^. The energy cutoff for the plane 
wave expansion of the PAWs was taken to be 550 eV. 
A (4x4x4) k-mesh has been used for self consistency. 
Symmetry has been switched off in order to minimize 
possible numerical errors. The calculation for the 
ferroelectric polarization with FM, AFM and NONC 
magnetic confmuration are carried out using Berry 
phase method [SI] as implemented in the Vienna ab initio 
simulation package (VASP). The Wannier function for 
the orbital (shown in Fig. 5) is constructed 



FIG. 6 . Different magnetic ordering (a) ferromagnetic (FM), 
(b) antiferromagnetic (AFM), (c) noncollinear (NONC). 


using the VASP2WANNIER and the WANNIER90 
codes. [ 1 ^ 


TABLE II. The relative energies and the change in Cu-Cu 
bond lengths upon relaxation for the nonmagnetic (NM) and 
different magnetic configurations with and without SOC have 
been listed here. -I- (-) signs indicate the increment (decre¬ 
ment) of the bond-length. 


Structure 

AE 

Change 

in bond lengths 


(meV) 

upon relaxation corresponding 




to the following 



exchange paths with respect to 



the experimental structure (A) 



J1 

J2 

J3 

J4 

J5 

NM-I-U relax 

130 

0.0 

- 0.01 

0.0 

0.0 

0.0 

FM-I-U relax 

120 

-0.07 

0.25 

- 0.01 

0.0 

0.0 

AFM+U relax 

82 

-0.06 

0.22 

- 0.01 

0.0 

0.0 

NONC-bU relax 

80 

-0.07 

0.23 

- 0.01 

0.0 

0.0 

FM+SOC-bU relax 

29 

-0.07 

0.25 

- 0.01 

0.0 

0.0 

AFM+SOC-bU relax 

0 

-0.07 

0.24 

- 0.01 

0.0 

0.0 

NONC+SOC-bU relax 

0 

-0.07 

0.23 

- 0.01 

0.0 

0.0 
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